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Abstract. The difference of vector and axial-vector charged current correlators is analyzed by means of 
QCD sum rules. The contribution of 10-dimensional 4-quark condensates is calculated and its value is 
estimated within the framework of factorization hypothesis. It is compared to the result, obtained from 
operator fit of Borel sum rules in the complex g 2 -plane, calculated from experimental data on hadronic 
r-decays. This fit gives accurate values of the light quark condensate and quark-gluon mixed condensate. 
The size of the high-order operators and the convergence of operator series are discussed. 
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1 Introduction 

The QCD sum rules PQ have been widely used for the de- 
termination of fundamental theoretical parameters, such 
as the coupling constant a s , quark masses and various 
nonperturbative condensates. Their accuracy depends on 
experimental errors and theoretical uncertainties. In many 
cases both experimental and theoretical errors are compa- 
rable by the order of magnitude, and any improvement is 
of interest. 

In this paper we will consider the 2-point correlators 
of charged vector and axial-vector currents, constructed 
from light u,d-quarks: 



= * J dxe**(TU}(x)U v (0)) 
= {q„qu-9„uq 2 )n^(q 2 ) + II$\q 2 ) 



where 



U = V,A 



The polarization functions iJW(s) have a cut along the 
real axes in the complex s = q 2 plane. Their imaginary 
parts (spectral functions) 



Vi/ffli(s) = 2irImII { yJ A (s + iO) 
a (s) = 27rImiT^ 0) (s + iO) 



(2) 



have been measured for < s < m 2 by ALEPH 2 and 
OPAL 3 collaborations from hadronic decays of r-lepton. 

Of particular interest is the difference Ily — iil , since 
it does not contain any perturbative contribution in the 
massless quark limit. The experimental data on the differ- 
ence vi(s) — ai(s) are shown in Fig.^ As demonstrated in 
0] , the dispersion relation can be written in the following 



form: 



i 

2V 2 



Vx(t) - ai(t) 



f 2 

dt + iz- 



E 

D>4 



O 



V-A 



\D/2 



(3) 



where the sum goes over even dimensions D of the oper- 

f 2 

ators (condensates) Od- The term /„• = 130.7 MeV 
is the pion decay constant, is the kinematical pole of the 
axial polarization function II A , see 0] for details. In © 

and below the notation 0^~ A stands for the condensates 
with all a s corrections, including slowly varying logarith- 
mic terms ~ ln™(— s). The list of the condensate contribu- 
tion to the vector and axial correlators separately can be 
found in 

The sum rules for the difference (01 have been studied 
in 0], 0-52 where the lowest order condensates 0% 
were found. Although the published values of 0^~ A are 
close to each other (within the errors), this is not the case 
for the operator 0^~ A . In 0], [Hj positive values of the 
D = 8 condensate were found, but the authors of recent 
publications [Hj, 52 have obtained negative condensate 
0\~ A . The source of this discrepancy could be very large 
condensates of dimension D — 10 and higher, accounted 
in ^21 : a typical ratio of the condensates in these pa- 
pers is |O2n+2/02„| ~ 5 — lOGeV 2 . If this statement is 
correct, the OPE analysis of 0] would be invalid, because 
the contribution of unknown high-order terms was esti- 
mated from the assumption \0\q/0§\ < 1 — 2GcV 4 . For 
this reason it would be interesting to find the operator 
O^Q~ A independently and compare it with the sum rule 
results. 

In this paper we repeat the analysis of 4 with the 
D — 10 operator included. In Section 2 all necessary V — A 
operators, obtained from the Operator Product Expansion 
in QCD, are listed and their values are estimated within 
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Fig. 2. Condensate expansion JHJ by the number of quarks in 
vacuum. Circles stand for the currents V/A, crosses are quarks 
in vacuum; gluons in vacuum are not shown. 



The contribution of the D = 8 4-quark condensates 
to the vector current correlator was originally obtained in 
|18| in factorized form and in |19| in complete (nonfactor- 
ized) form. In 0] these results were verified and an am- 
biguity of the factorization at the N~ 2 order was pointed 
out. Here we will follow the factorization procedure, de- 
scribed in Appendix B. The result is 1 : 



O 



V-A 



= ?>-KC N a s ml(qqY 



(6) 



where the mass mo is defined from the 5-dimcnsional quark- 
gluon mixed condensate: 



the framework of the factorization hypothesis. In Section 3 
the operator values are obtained from the fit to Borel sum 
rules. In the last Section the validity of our assumptions 
is discussed and the results are compared with the ones 
obtained in other publications. The complete form of the 
D = 10 operator and technical details of its derivation are 
dropped to Appendices A,B. 



2 V-A operator expansion 

The first term in the operator series is the D = 4 
operator: 



0\~ A = 2(m u + m d )(qq) 



1 + 4 a s (Q 2 ) | 59 / a s (Q 2 ) 

3 7T 6 \ 7T 



(4) 



where Q 2 = ~q 2 , we assume (uu) = (dd) = (qq). The a s 
corrections have been computed in U2| an d UU- in iac t, 
the contribution of the D = 4 operator to the sum rules 
considered here, is small. So we can safely neglect the a s - 
corrections in (@J and put C\~ A = -f 2 m 2 = -3.3 x 
10~ 4 GeV 4 , as follows fr om Gell-Mann-Oakes-Renner low 
energy theorem |15| . 

The D = 6 operator in factorized form is equal to: 



a 



V-A 



= -8wC N a s (qqy 



1 + 



Ap 2 ) 



C6 



7T \ 4/1^ 

(5) 

where Cn = 1 — N c 2 = 8/9 is the color factor, which 
appears in the factorization of the 4-quark operators at 
the leading a s -order. The NLO terms were computed in 
[If)] and the constant cq was found equal to 247/48. In UJj 
another treatment of 7 s matrix in dimensional regulariza- 
tion was employed, leading to eg = 89/48. For the later 
choice at fi — 1 GeV and a s (/i 2 ) = 0.5 one finds the fac- 
tor in square brackets in (JjJJ equal to 1.3 (the logarithmic 
term can be neglected due to small numerical coefficient). 



l(qGq) = 2(qD 2 q) 



■m%(qq) 



(J) 



where G = jaJpGap, G a p = i[D a ,Dp] is the gluon field 
strength, sec Appendix A for more definitions. The pa- 
rameter ttjq has a meaning of typical momentum of vir- 
tual quarks in vacuum. It was found from baryonic sum 
rules ml = 0.8 ± 0.2 GeV 2 |2oT2T] . and B - B* splitting 
[2"2] , The values close to 1 GeV 2 were also obtained from 
the latest lattice calculation 23] and in QCD string model 

pa- 

There are many different condensates of dimension 
D = 10. They can be grouped into four parts: 



nV-A _ q 



(0) 
10 



O 



(2) 



o 



(4) 



o 



(6) 



(8) 



where upper index (i) denotes the number of quarks in 
vacuum. This separation is shown diagrammatically in 
Fig. |21 The purely gluonic operators and the 2-quark 
ones cancel in the V — A correlator in the limit of 
massless it, <i-quarks. The operators with 6 quarks in vac- 
uum have the structure ((qq) 2 (q~Dq)) . After factorization 
they become ~ m(qq) 3 , which is again negligible for light 
quarks. The only essential contribution to the V — A sum 

rules comes from the 4-quark operators . 

In this paper we have computed the contribution of 
the 4-quark condensates to the vector and axial current 
correlator. Details of calculation and complete form of the 
operator are given in Appendix A. The factorization 
scheme necessary to reduce large number of independent 
structures, is described in Appendix B. The result is: 



O 



V-A 
10 



N 



T <<?Gg> 2 



16 ( 3Xi - X 2 + X 3 + -Xi ) (qq) 



(9) 



1 In |1| the factor Cjv was ignored, since 0(N C 2 ) terms were 
neglected 
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where Xi are 4 independent D 
sates: 



7 quark-gluon conden- 



Xi = (qGapGapq), 
X3 = {q r YapG a - f Gi3~ / q), 



X 2 = i{qi b G a pG al 3q), 
X4 = i(qj a p{D a Jp)q). (10) 



where G a p is the dual gluon- field strength, 7^/3 = 5(7,37/3— 
Ipla) and J a = DpGap. Their numerical values arc not 
known. The condensate X4 can be brought to the 4-quark 
form X4 ~ ((qq)(qDq)) ~ m(qq) 2 , which is negligible. In 
order to estimate other condensates, we assume further 
factorization according to (qTq) — (qq) (tr _T)/(4iV c ), the 
trace is taken both over color and spinor indices. Then: 

X x = — a s {G a afj G a al3 ) (qq) , X 2 = X 3 = X 4 = . 

(11) 

Under these assumptions the operator takes the form: 

50 



O 



V-A 



■ira s C N (qqY 



■ ml + 327ra s <G 2 ) 



(12) 



It is rather difficult to find accurate value of the gluon 
condensate from any sum rule. Detailed analysis of char- 
monium sum rules performed in |25| has lead to the re- 
striction (^G 2 ) = 0.009 ± 0.007 GeV 4 , in agreement with 
many previous estimations. Taking this central value and 



1 GeV^, one obtains 0\^ A jO 



V-A 



\V-A 
J 6 



0.8 GeV . For 

$±2.1) x 10~ 3 GeV 6 0] we find the following 



estimation of the D = 10 V — A condensate: 



O 



V-A 
10 



= -5 x 10~ 3 GeV 



10 



(13) 



In the next section we will compare this estimation with 
results of the fit, obtained from the sum rules. 



3 V-A sum rules 

Many different sum rules have been investigated in order 
to determine numerical values of the condensates. Most 
of authors employ polynomial sum rules: the correlator 

n$\s) - II%\s) is multiplied on some polynomial of s 
and then integrated over the circle \s\ — sq in the complex 
s-plane. Their advantages are: 1) one does not need to 
know the spectral function vi(s) — a±(s) for s > so, which 



allows to reduce high error from the region s 



by 



choosing so reasonably below m 2 and 2) all operators of 
dimension higher then the polynomial dimension, do not 
enter these sum rules due to Cauchy theorem. But the 
disadvantages are also obvious. If the operator expansion 
(j2J is divergent (asymptotic), the Cauchy theorem is not 
applicable to this series. Moreover, possible logarithmical 
terms ~ ln k Q 2 /Q 2n appear at the NLO in a s expansion. 
These terms contribute to any polynomial sum rules. It 
makes uncontrollable the contribution of the high order 
operators to the polynomial sum rules at so < 2 GeV 2 , 
especially for large ones as obtained in [9TT2"] . 

For these reasons we prefer Borel sum rules, where the 
high order operators are suppressed as Oi n ln\. In order to 



separate out the contributions of different operators from 
each other, one may consider the Borel transformation in 
the complex plane of the Borel mass M 2 — ► M 2 e l ^~^ 
(which is equivalent to the Borel operator applied to the 
dispersion relation @ written along the ray s — > se 1 ^ in 
the complex s-plane 0]). The real and imaginary parts of 
the Borel transformation are: 



= / 2 + ^(-) fcC ° S(fc ^ 02/fc+2 



k=l 



kl M 2k 



exp 



(m 2 



cos<^ sin y-j^ smc^j (vi - ai)(s) 



-Eh 



. (k4>) o 



V-A 
2k+2 



k=l 



ft! M 2k + 2 



(14) 

ds 
2n 2 M 2 

(15) 



We made the imaginary part 1)15(1 dimensionless, while the 
real part (|14fl has dimension GeV 2 in order to separate 
out the leading constant term f 2 . The logarithmical terms 
are neglected in the rhs of I(14I15JI . otherwise the terms ~ 
In M 2 appear. The only known logarithmical term is in the 
a s -correction to the D = 6 operator It can be easily 
taken into account (see [5] for explicit formulae), but its 
relative contribution is negligible due to small numerical 
factor, so we shall ignore it. 

The derivation of (|14I15|1 from the dispersion relation 
Pjl implies infinite upper integration limit s m = 00. Ex- 
perimental data on the axial function ai(s) are available 
only for s < m 2 = 3.16 GeV 2 . However the data at s > 
3 GeV are rather unstable and have large error because 
of low statistics, see Fig ^ For this reason we put s m = 
3.0 GeV 2 in H14I15[) . Removal of the data above this point 
does not change the Borel transform significantly (if M 2 is 
not sufficiently large), but may reduce the errors. In fact, 
the sum rules considered here do not rely on the high- 
energy data: say, if the upper integration limit s m is re- 
duced to 2.5 GeV 2 , the condesates change at most within 
10% limit. If the data above 3 GeV 2 are removed, both 
ALEPH and OPAL data give almost equal central values 
and similar errors of the Borel transforms (|14I15|1 . For this 
reason we will psesent below the analisys of ALEPH data 
only, since they have smaller errors. The condensates, ob- 
tained from OPAL data are almost the same. 

The argument of the exponent must be negative cos 4> < 
in order to suppress contribution of the high-energy 
states from unknown region s > m 2 , which means ir/2 < 
cj> < 7r. Of special interest are the closest to 7r (minimal 
error) angles at which the contribution of some opera- 
tor 02/c+2 vanishes. Such angles are (f> = n(2k — l)/(2fc), 
k = 2, 3, . . . for real part (|14J) and <fi — ^(k — l)/ft, k = 
3,4,... for imaginary one l|15fl . The sum rules (|14I15H 
at some of these angles were considered in 0] with the 
operators Oq and 0$ as free parameters to fit. It was 
shown, that for O v ^ A = -(6.8 ± 2.1) x 10" 3 GeV 6 and 



O 



V-A 



(7 ± 4) x 10~ 3 GeV they are well satisfied for 



M 2 > 0.6 GeV 2 
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M 2 , GeV 2 


n v ~ A n v ~ A n v ~ A 


Xo 


27r/3 


0.4-1.0 
0.5 - 1.0 


-7.0 ±1.4 -3.8 ±3.3 
-7.0 ±2.5 -3.9 ±10.3 


0.14 
0.14 


3tt/4 


0.4-1.0 
0.5 - 1.0 


-7.3 ±1.1 8.0 ±2.1 
-7.9 ±2.3 9.5 ±2.5 


0.17 
0.12 


47r/5 


0.3-1.0 
0.4-1.0 


-8.1 ±1.8 10.3 ±4.1 -7.6 ±5.0 
-9.3 ±5.1 13.8 ±15.1 -13.2 ±23.8 


0.11 
0.06 


all 




-7.2 ±1.2 7.8 ±2.5 -4.4 ± 2.8 
-7.5 ±2.3 8.6 ±6.0 -5.2 ± 8.4 


0.40 
0.20 



Table 1. Operator fit obtained from eq 115H for different angles <f>; operators are in 10 3 GeV D . Last two lines contain combined 
fits for all these angles for upper/lower choice of M 2 range. 



It is more difficult to find high-order condensates (say 
0\i A ) from the sum rules, since several unknown param- 
eters enter the same equation and the high-order conden- 
sate strongly depends on exact values of the low-order 
ones. Here one needs to consider the Borcl transformation 
at several values of M 2 , where the relative contributions of 
various condensates are different. In other words, we may 
fit the shape of theoretical curve with experimental one 
within some reasonable region M 2 < M 2 < A/f • For this 
purpose it is natural to define the least square deviation, 
normalized to experimental error: 

where B th eor / B ex P is the right/left hand side of the Borel 
sum rules (I14I15|I . One may calculate x 2 with theoretical 
condensates Oi as free parameters. It is quadratic function 
of them: 

X 2 = xl + E C « (Pi - i°i ~ 3 ) (17) 

Obviously Oi are the central values of the condensates. 
According to the definition Ijlfijl it is natural to consider 
equation x 2 = \ as the one, which determines the border 
of the ltr deviation area in the parameter space. Diago- 
nalizing the matrix Cy by means of orthogonal rotation 
we conclude, that Cjj is inverse to the covariance matrix 
AOi ■ AO. } = ((7-%-. For a good fit xl < 1- 

The fit results depend on the Borel mass limits M 2 2 

in (|16fl . For M 2 > 1 GeV 2 the experimental errors are 
large, so we take M 2 — lGeV 2 . The lower limit M 2 
depends on the size of neglected high order operators. 
In 0] a good coincidence of experimental and theoreti- 
cal curves was observed for M 2 > 0.6 GeV 2 . Here we in- 
clude the operator O\o in the analysis, so this value can 
be slightly reduced. As follows from our calculation of the 
4-quark condensates J5J, 10, 1|12|) . it is reasonable to as- 
sume C>2n+2/02n ~ tuq ~ 0.7 GeV 2 . It leads to an es- 
timation |Oj^T A | - 3 x 10~ 3 GeV~ 12 which allows us to 
take Mi = 0.4 GeV 2 , where typical contribution of such 
operator is not higher than 20%. At the angles, where the 
contribution of the operator On vanishes, the Borel mass 
can be reduced even further, say, to M\ = 0.3 GeV 2 . All 
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Fig. 3. Confidence level contours, obtained from 2-parameter 
fits of the sum rule l|15^ in the range M 2 = 0.4 - lGeV 2 . 
The condensates are in 10 -3 GeV D , the contours show 1, 2, 3<r 
deviations of \ 2 ■ 



these assumptions are confirmed by the results of the fit, 
see Figures below. 

The condensates, obtained from real part of Borel trans- 
formation (|14fl . are sensitive to exact value of /„•. For this 
reason we shall use the imaginary part (|15Jl for numer- 
ical fit. The best angles are 4> = 27r/3, 37r/4, 47r/5 where 
the contribution of the operators Os, Oio, O12 vanishes re- 
spectively. The fit results for each angle are summarized 
in the Tabled The lowest errors are obtained from the 2- 
parameter fits at the first two angles. The deviation Xo for 
these fits is sufficiently small. For this reason the inclusion 
of additional parameters, say O12, will not improve the fit 
quality, but will increase the errors only. 

The operator values, obtained from the sum rules, are 
not independent but have large covariances 

Pii = AO~AD~/((AO~j 2 jAOtf) 1 ' 2 . 

All fits give p6,io ~ f and pe,a ~ Ps.io ~ — 1- For the 2- 
parameter fits the covariances can be demonstrated on the 
confidence level plots, see Fig EJ The equations x 2 = n 2 
set the ellipses, which are the borders of the ner deviation 
area. 

One may also try to fit the condensates at all these 
angles simultaneously by minimizing x\u = |[x 2 (2tt/3) ± 
X 2 (37t/4) ± x 2 (47r/5)], see the last two lines in the Table. 
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Fig. 4. Imaginary part of the Borel transformation for <f> — 
2n/3 (no Og) and (j> = 37r/4 (no Oio). Shaded area is the 
l.h.s. of 1151 calculated from experimental data (with error). 
The lines display the operator series in the r.h.s. of I15H with 
condensates, equal to the central values of l|18jl . The number 
nearby each line shows the order of the series; say "8" denotes 
the contribution Oi+Og+Os. Grid shows possible contribution 
of the operator 12 within the limits \0^f A \ < 3x 1CT 3 GeV 12 . 



As the final result of our analysis we take this combined 
fit: 



V-A 
6 

V-A 
8 

V-A 
10 



- (7.2 ± 1.2) x 10" 3 GeV' 
(7.8 ±2.5) x 10~ 3 GeV 
-(4.4 ±2.8) x 10" 3 GeV 



■6 
•8 
10 



(18) 



The lower limit of the Borel mass in (|16f) was taken M 2 — 
0.4 GeV 2 for the first two angles and M? = 0.3 GeV 2 for 
the last one. If M 2 is taken by 0.1 GeV 2 higher, the errors 
are increased, especially for the high dimension operators, 
see the last line in the Table. 

The validity of our assumptions is demonstrated in the 
Figure 0] If the operator 0\^~ A is taken into account, a 
good agreement of theoretical and experimental values is 
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Fig. 5. Imaginary part of the Borel transformation l)15|l at 
4> = 4tv/5 and real part 114H at <j> = 97r/10. The operator 
0X2 ~ A vanishes in these sum rules. 



observed for M 2 > 0.4 GeV 2 . Below this value the con- 
tribution of the operator Oyi could be large. Even better 
agreement can be found at the angles, where the oper- 
ator O12 disappears, see the plots in the Figure |3J Here 
the fit can be extended down to M 2 — 0.3 GeV 2 . One 
may also obtain the condensates by fitting the real part 
of the Borel transformation l|14|) . Here the central values 
of the condensates turns out to be close to l|18(l . but the 
errors are higher due to the presence of additional pa- 
rameter / 2 . Combined fit of eq l|14|) at different angles <fi 
gives = 131 ± 4MeV. As pointed in 0], itself has 
an ambiguity of order m 2 /m 2 ~ 3%, the accuracy of the 
chiral lagrangian parameters. Notice the sign alternation 
in (|18(l . in agreement with the minimal hadronic ansatz 
for Ily — 11a correlator, constructed in |7| in the large N c 
limit. 

Finally, we write down the values of the quark conden- 
sate and the parameter m§, obtained from the operators 
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GE}: 



a s [qq)- [l + c e — ) = (262 ±9MeV) t 



(19) 



= -ol 



l /ol 



= l.l±0.2GcV 2 (20) 



The errors in the r.h.s. are purely experimental: they do 
not include possible contribution of the operator On and 
higher as well as unknown QCD corrections to the con- 
densates. The factor c§ is scheme dependent and left arbi- 
trary in l|19|) ■ The accuracy of m§ is better than accuracy 
of 0^~ A because of high covariance of Oq and 0%. Notice 
a very good agreement of the D — 10 condensate (|18fl 
obtained from the sum rules, with the one estimated in 
previous section i|13|) within the framework of the factor- 
ization hypothesis. 



4 Conclusion 

We have performed the analysis of the V—A spectral func- 
tions, obtained from hadronic r-decay channels, with the 
help of the Borel sum rules. The values of the condensates 
of dimension D = 6,8, 10 were found l|18(l by fitting the 
theoretical curves of the Borel transform to the experimen- 
tal ones within its errorbands. The major contribution to 
these condensates comes from the 4-quark operators. Its 
contribution to the current correlators was calculated and 
their size was estimated by means of the factorizations hy- 
pothesis. The estimated value of the D = 10 condensate 
(|13|) is found to be in good agreement with the fit result 
<|18[) , which demonstrates the validity of OPE approach in 
Quantum Chromodynamics. 

Our results are based on several assumptions, in par- 
ticular, the factorization (vacuum insertion) hypothesis. 
There is a statement in literature |26) . that factorization 
hypothesis underestimates the quartic condensates by a 
factor ~ 3. This conclusion is based on the comparison 
of the quark condensate obtained from the D = 6 opera- 
tor in yO-meson (vector) sum rules with the one calculated 
from the low-energy GMOR theorem. (Our result (|19fl is 
also larger than GMOR condensate for reasonable theoret- 
ical parameters.) However this comparison has many other 
sources of error, such as scale-scheme ambiguity, high or- 
der QCD corrections, light quark masses, corrections from 
the chiral lagrangian etc. The accuracy of the factoriza- 
tion hypothesis can be of the same order as the ambigu- 
ity of the factorization of the D=8 operators the level of 
0(N~ 2 ) terms, as demonstrated in 0]. More careful state- 
ment about validity of the factorization hypothesis could 
be obtained by evaluating the contribution of the meson 
states to the 4-quark condensates. 

Second objection may concern rather low value of the 
Borel mass Ml = 0.4 GeV 2 used in our fit (|16|) . Indeed, 
the typical scale where perturbative results for the current 
correlators are confirmed, is Q 2 > 1 GeV 2 . But our result 
for the D = 10 operator demonstrates rather low (power- 
like) growth of the operators |02n+2/02 n | ~ n^o in the 
V — A channel. If the operators grow as \02n\ ~ tUq, then 



the Borel series behaves as e m o/ M . The contribution of 
the n+l-term in the exponent is small for M 2 ^> m 2 /n. So 
for n = 5 the minimal scale M 2 = 0.4 GeV 2 seems reason- 
able. For a faster growth of the operator series this choice 
could be inappropriate. For instance, if one plots the Borel 
transformation versus M 2 with the condensate values ob- 
tained in the divergence of the operator series will be 
obvious already at M 2 w 0.7 GeV 2 . However it should be 
mentioned, that the D — 10 condensate obtained there 
exceeds our value (|13|) by an order of magnitude. It seems 
unlikely to explain such discrepancy by the inaccuracy of 
the factorization. All these assumptions can be confirmed 
or disproved only within a nonperturbative approach. 

We have neglected the logarithmic terms ~ hiQ 2 /Q 2n 
in the OPE series ©. Such contribution from the as- 
correction to the D = 6 condensate © has a small nu- 
merical factor; its discontinuity along the real axis Q 2 = 
— s < is too small to compare with the spectral function 
Vi(s) — oi(s). For this reason it would be interesting to 
calculate the a 2 correction to the D = 6 V — A operator 



and a s correction to the operator 8 
in the sum rule analysis. 



V — A 



and include them 
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Appendix A: 4-quark operators 

The calculation of the operator contribution to various 
current correlators can be performed within the frame- 
work of background field method, see for instance |27| . 
Here we describe the algorithm, conventions and basic 
formulae, necessary to calculate the contribution of the 
4-quark condensates to the 2-current correlator, which 
correspond to the third diagram of the Fig. |3 We also 
present here complete form of the 4-quark operators up 
to dimension D = 10. For definiteness we consider only 
the vector current correlator; the condensate contribution 
to the axial current correlator is trivially obtained by the 
substitution d — > j 5 d. The contribution of the 4-quark 
condensates can be written as: 



nUq) 



' 4 



dx dy dz e 1 



D ab 



(Al) 



uix^Six, y)-y a X a d(y) + u(y)^ a X a S(y, x)-y^d(x) 
d(0) 7 „S(0, z) ll3 X b u(z) + d(z) ll3 \ b S(z, 0) 7 „u(0) 



Here S(x,y) — (Tq(x)q(y)) is the quark Green function 
and D^(x,y) — (Ta'J l (x)a b J (y)) is the gluon Green func- 



tion in background gluon field A^ 
the equations: 



A^+Qfj,. They obey 



iD x S(x,y) = iS 4 (x-y) (A2) 

[D 2 x9m + 2G m {x)} ac Dt(x,y) - i5 ab g ^ (x - y) 

(A3) 
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where g^ v = (+,—,—,—) is Minkowski metric. The quarks 
are massless, the gluon Green function is taken in the 
Feynman gauge. The covariant derivative and the gluon 
field-strength tensor in the fundamental representation 
<|A2(1 are defined as follows: 



Gfj,v = i[DfM,D v ] = -X a G^ v 



where A a are Gell-Mann matrices tr{\ a \ b ) = 2S ab , [A a , X b ] 
= if abc \ c . We shall also use additional compact notations 
for these objects in adjoint representation (IA3I) : 

Df = d ll 5 ab +A«», Af = 9 -f acb Al, Gf u = p acb G^ 



It is convenient to perform partial Fourier transforma- 
tion of the Green functions: 



S{x,y) 



d 4 q 



-iq{x-y) 



(2tt)4 

d 1 -iq(x 
(27T) 



S(q,y), 



e- iq{x - y) Df v {q,y) (A6) 



Then one can write down the solution of equations (|A2I 
(j A3|l as series in powers of background field A: 



oc 



n=0 



(A7) 



DfMv) = {a>(«)£ [iR(q,z)D (q)} n }?», (as) 



?i=0 



where 



So(q) = -o: 



D (q) 



are free propagators, x = y — id 7 the derivative d — d/dq 
acts on everything from the right as [ d ^, q v ] = <7 M „; R is 
the following matrix operator: 

= [~^ a A ab {±) - iA ab {x)q a 

+A a a c (x)A c a b (x)]g fW + 2G; b / (x) (A9) 

The equations (|A7|) . I|A8() can be evaluated in a gauge 
covariant way in the fixed point gauge x^A^{x) = 0, where 



da olx v G^yiax) 



(A10) 



E 



(n + 2)n! 



D ai ...D an G^(0) 



q{x) = £ : 

n=0 



D ai ...D an q(0) 



(All) 



In order to compute the propagators S, D for any fixed 
order n, one has to substitute (|A10|) into (|A7(I . l|A8Jl . move 



all the derivatives d to the right and then leave only the 
terms without d . 

The 4-quark condensate contribution lAlfl can be writ- 
ten in terms of the propagators S, D as follows: 



(A4) nlM 



u{~i d hf,S(q, -i d )~f a \ a d(-i d ) 



x D%(q, -i d ) X b (q) + X b fj (q) D b f a (-q, -i d ) 



x u(-i d )-f a X a S(-q, -i d h„d(-i d ) 



(A12) 



(A5) 

where 



Y b 



(, l) =d{-i~d) lp \ b S{q, 0) 7 ^(0) 

f d(0)7„S(-?, -it)j \ b u(-i*d) . (A13) 



In the functions S(q, y) and D(q, y) the derivatives d , d 
over momentum q always stand on the right from any func- 
tion of q: ... q ... d. The derivatives inside (|A13|I do not act 
on anything outside X b ^ . After these derivatives are evalu- 
ated, we compute the transverse part 77^ = —Il^^/iZq 2 ) 
defined according to Q). (We also checked, that longitu- 
dinal part vanishes 11^ = 0.) And finally, to separate out 
the Lorentz invariant condensates, we average over 
directions of vector q^ according to: 



(2n-l)" 2 



■ <W = 2 



QfJ.1 ■ ■ ■ q^2n+l 



(2n + 2)!! 
0, 



(A14) 



where (fii . . . /i„) denotes usual index symmetrization with 
weight 1/n! . All these calculations were performed on com- 
puter. 

The most time-consuming part of the calculation is to 
reduce large number of terms in the final result to min- 
imal number of independent structures. For this purpose 
we employ the quark equation of motion Du = Dd = 
and the "integration by part" identity (A(D^B)) = 
~{(D^A)B) (the vacuum average of the total derivative is 
zero (d x O(x)) — d x (0(x)) = d x {O(0)) — for any gauge 
invariant operator 0{x)). It allows to bring the operators 
to obviously hermitean (real) form, which provides an ad- 
ditional verification of the result. 

In order to write down the 4-quark condensates in a 
compact form, we introduce here the following bilinear 
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quark structures of increasing dimension D: 

3D: A a = d~X~f 5 ~f a u 

AD : B^p = i(d\j a up - dpX-y a u) 

5D : = i{d^l»Uf3j ~ d^XjaU) 

C al-t = rfA7 5 7 Q U/3 7 + d ( 3 7 A7 5 7 Q u 
6D : E^l s = d{\, G Q/3 }7 7 U5 + ds{X, G Q/3 /7 7 u 

(2) — ~ — ~ 

E af3~/8 = Gapjj^Us + d S {X, G al3 }jjU 

Ea/SjS = i(d{X, G Q/3 }7 5 7 7 u 5 - d 5 {X, G a/3 }7 5 7 7 u) 



'a/3-yS = i{ d {\ G q/3 }7 5 7 7 m 5 - d s {\, G Q/3 }7 5 7 7 u) 



E a0-yS = rf{A,G /37;Q }7 (5 U 
7D 1 F a/3-y = ^a}7 5 7/3M 7 ~ d ^{X, J a }l b lf]U) 

*Sr = ^ i G ^, Gs^W^u 

F a/3~/ = id {^i i G aS, G < 5 /3 ]}7 7 U 

(A15) 

where u a — D a u, u a p = D^D^u = \{D a D p + D pD a )u, 
G/3 7;Q = DaGpj, J a — DpG a p, [A, B] = AB - BA, 
{A, B} = AB + BA. The dual tensor is defined as G a p = 
^EaPuvGfiv, e 0123 = 1 and 7 5 = i^^^j 3 . The values 
G, G, J in (|A15p are in fundamental representation. All 
bilinear structures belong to adjoint representation of the 
gauge group, the gauge index of the Gell-Mann matrices 
A is omitted. We denote conjugated structures by over- 
lined letters, which are simply obtained by the replace- 
ment u <^ d, for instance A a = A^ a = uXj 5 j a d. 

The 4-quark condensates of dimension D = 6, 8, 10 are: 



0% = -2ira s (A a A a ) 



(A16) 



01 



A d(1)d(1) r( 2 )r( 2 ) 



-ACf a(i A a - AAJjf^ + l2A a G a pAp) (A17) 



OX Q = H±( 25G« ,cil. - 5&±C*L - iog<1g?> 
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a/?7 a/?7 afi-y a/37 



J (xficx. 7/37 



lar^ 4 ) r( 4 ) _ 15 r<( 4 ) _ ar^ r<- 4) 

- 2B ap G Pl B a] ~ 66B^G aJ B^ + A a (8J[ a] 0] 

\ 33 - 

— 3G Q7 G 7/ 5 + 19G / 3 7 G 7Q J Ap + —A a Gj3 1 Gi3 1 A a 

, o(l) , 5 p(4) _ 7 P (5) _ oor< o(2) 

' a/3 I /9707 2 "77/3 2 7' 3 7 Q ^^Pl^ory 

4- 21 r< 4 V f pW _l 5 f< 4 ) 7 p( 5 ) 

T ^ ^^7; "^7 I ~r I £/ /3 7 a 7 T 2 ^ajyfj 2 lPl a 

+2&B<$G fh - yI T G ft;a j flg + fig (-H^C 



(2) 

77/3 



15^(2) ^ 5 ^(3) 



(3) 



1 



(6) 



(6) 



-4J5. 



(6) 



17 



7 a/3 7 2 "°>w 7^ 



G Q/ 3 ;7 A 7 — llJa^Q, 



11 



(2) 



2 ^0:77/3 



4 0ia-f 2 "7/37 J& /37a7 "T" 2 0/377 2 /3ce7 ' 7 



-4K 



(6) 



17 



7a/37 2 



+ HA a J/3 73^' 



?(2) 



o I or 3«3 + W« + 2 /3/3q ( a /3)/3 



1 



(2) , -(3) 



_op(1) , , 9 p(2) _ 1 p(2) , p(3) 



, 45 ^(4) ^ 35 - (4 ) X 75 - (4 ) , 

+ Y C/ /377 Ct "/3 ~ T oft ^ ~ ~°/37Q C:r /37 ) A " , 



2 a/37 g /37c 



(Ais) 



In (|AT7|) . (|AT8|) the field strengths G, G, J are in adjoint 
representation G™« = | f acb G c a p etc; gauge indices are 
omitted, say A Q G/3 7 G / g 7 A Q denotes A^Gf^G^A^. The 
operator Og" (|A17|I can be easily brought to the form, 
obtained in 0j and 



Appendix B: Factorization of 4-quark conden- 
sates 

At first let us remind, how the factorization (vacuum in- 
sertion) works for the D = 6 operators. It is illustrated by 
the following equation: 

((u\r 1 d)(d\r 2 u)) = -2Gjvtr [(u(g>u}ri(d<g>J)r 2 ] 

(Bl) 

where /Ti are some Dirac matrices, G/v = 1 — l/A^ 2 , iV c 
is the color number, kept arbitrary here. In (IB 1|) the no- 
tation (q (£) q) denotes 4x4 matrix in spinor space, the 
color indices are contracted. It is proportional to the quark 
condensate: 

(q ® <l) = - t («9> (B2) 



The result of the factorization is well known: 



01 = ~4:Tra s C N {uu)(dd) 



(B3) 



(In the vector sum rules one also accounts additional op- 
erator (q"f a J a q), which takes the 4-quark form when the 
gluon equation of motion is applied. Such operator comes 
from the 2-quark diagram, so it cancels in the V — A cor- 
relators.) 

The factorization procedure becomes ambiguous at the 
level of D = 8 4-quark condensates. As shown in fj], dif- 
ferent ways of factorization give different terms ~ V-^c ■ 
For definiteness, here we follow the following factorization 
scheme. At first we replace the field strength by the deriva- 
tives as G^v = ilD^D,,] for fundamental representation 
and G ab u = D„] ab for adjoint one. Then we apply the 
equation I|B1J1 . where the quark wave functions u, u and 
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d, d may carry some derivatives. Finally, the quark ma- 
trices (...) with derivatives are expressed in terms of the 
condensates as: 

(D a q ® q) = 0, 

{D a Dpq ®q) = ~ (a a p + ^afi) i(qGq), (B4) 

where 7 a/3 = 7^7,3] = 5(7^ 7/3 ~ Ipla), G = j a pG a p. The 
result for the D = 8 condensate is: 



-27ra s C 



AT 



(uu)i(dGd) + (dd)i(uGu) (B5) 



In the condensate Oio one encounters the terms with 
quarks carrying 4 derivatives. We average these terms with 
the help of the following rule: 

(D a D p D 1 D s q®q) = -—[g a ^ s (3X 2 +6X 3 -2X 4 ) 

+9a- l 9l36{6X 1 + 3X 2 + 6X 3 + 4X±) + g aS g Pl {12Xi 
+3X 2 + 6X 3 + 4X 4 ) + {g a p llS + g lS -f a p)(X 1 + 2X 2 
+3X 3 ) + (ga-yjps + g sj ai )(2X 1 +X 2 + 3X 3 + X 4 ) 
+g a sjp 1 {X 1 + 2X 2 + X 4 ) + g fj ^ aS (Xi + X 4 ) 
+3j af}lS X 2 ] (B6) 

where jap-yS = 7[a7/3777(5] j Xi are 7-dimensional conden- 
sates, defined in ijTU|) . 

Being applied to the operator l|A18|) . this procedure 
gives the following result: 
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(uGu)(dGd) 



7 



-4(3Xr-X 2 « + X 3 " + -X 4 ") (dd) 



4 | 3Xf - X$ 



xs 
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where Xf are constructed from the quark of flavor q. Fhe 
axial condensates can be obtained by simple replacement 
d "f 5 d. For all factorized 4-quark operators O^, = — O^. 
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